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' Abstract 



In this note, we present a new numerical method for solving backward stochastic differential 
equations. Our method can be viewed as an analogue of the classical finite element method 
solving deterministic partial differential equations. 



Qh ! 1 Introduction 

sz . 

Linear and nonlinear Backward Stochastic Differential Equations (BSDEs in short) were introduced 
in [1] and [9], respectively. It is well-known that BSDE plays crucial roles in Stochastic Control, 
Mathematical Finance etc. Clearly, for applications, it deserves to develop effective numerical 
methods for BSDEs. 

\ Let T > and (0, J 7 , F, P) be a complete filtered probability space with F = {Ft }te[o,T] j on which 

a 1-dimensional standard Brownian motion {w(t)}tg[o,T] 1S 

defined. We denote by L 2 jrt (0;M. n ) (n G 
N) the Hilbert space consisting of all J-^-measurable (R n -valued) square integrable random variables; 
by L|(O;L r (0,T;lR n )) (1 < r < oo) the Banach space consisting of all M n -valued {J^j-adapted 
| processes X(-) such that E|X(-)|^ T . R n) < oo; and by £f(0; D([0, T]; W 1 )) the Banach space 

consisting of all M n -valued {J^j-adapted cadlag processes X(-) such that E(|X(-)|? 



L|°(0,T;R™)- 



< OO. 



For y T G L 2 TT (n-,M. n ) and /(-,-,•) satisfies /(-,0,0) G L 2 (J7; L 1 (0, T; M n )) and the usual globally 
Lipschitz condition, we consider the following BSDE 



>< ■ ( dy(t) = f(t,y(t),Y(t))dt + Y(t)dw(t) in [0,T], 



d ■ I y(T) = yr 



(1.1) 



Various numerical methods have been developed to solve equation (jl.ip . say in [21 [TJ [TUl H] and the 
references therein. These methods use essentially the strong form of (|l.lj) . which holds true only 
if F is the natural filtration generated by the Brownian motion. Also, it seems that the previous 
methods need to compute the conditional expectation, which is in general not easy to be furnished 
numerically. 

In this Note, we shall present a new numerical method solving BSDEs from the viewpoint of 
transposition solution introduced in [6], as recalled below. 
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Definition 1.1 A couple (y(-),Y(-)) € L|(0; D([0, T];R n )) x L|(fi; L 2 (0, T; R n )) is caZZed a irans- 
position solution of BSDE {73}, i//or any f € [0,T], («(■),«(■), 77) G L|(0; L x (i, T; M n )) x L 2 (Q; 
L 2 (i, T; M n )) x L?^ (f2; R n ), i/ie following variational equation holds 

E(z(T),y T ) -E( V ,y(t)) 

= E f (z(r),f(r,y(r),Y(r))dT + E f (u(t), y{ T ))dr + E / <„(t), F(r))dr, (1.2) 

Ji Jt 

where z(r) = rj + J t u(s)ds + J" 4 T v(s)dw(s). 

We refer to [6] for the well-posedness of equation (jl.ip in the sense of transposition solution. It 
is easy to see that, if this equation admits a strong solution (say under the assumption of natural 
filtration), it coincides with the transposition solution. 

Our numerical schemes for solving equation (II. ip can be described as follows. 

1) Take a suitable finite dimensional subspace H m of L|(r2; L 2 (0, T; M n )); 

2) If the solution of (jl.ip exists, then it should satisfy the variational equation (jl.2p for any 
u, v £ -ff m . By taking it and v to be the orthonormal basis of H m , we obtain a system of 
approximating equations; 

3) If the solution of the system of approximating equations exists, then we find a class of nu- 
merical solutions of (|1.1D : 

4) Finally, we show the convergence of the above numerical solutions. 

Clearly, the above procedure is, in spirit, very close to that of the classical finite element method 
solving deterministic PDEs (e.g., |3j). Therefore, our method to solve BSDEs can be viewed as a 
stochastic version of the finite element-type method. Nevertheless, the notion of "stochastic finite 
element method" has already been used for other purpose, say [U [8] and references therein for 
solving random PDEs. Note also that our method is quite different from that in these references, 
and therefore instead we call it a finite transposition method. 

There are at least two reasons for us to develop this new numerical approach for BSDE. The 
first one is that, we can solve the BSDE with general filtration. The second is that, in our approach, 
we do not need to compute the conditional expectation. 

We refer to [11] for the details of the proofs of the results in this Note and other results in this 
context. 



2 Numerical schemes and convergence 

For simplicity, we consider only the following linear BSDE (with /(•) £ L^.(0,;L (0, T;R n ))) 

/ dy(t) = f(t)dt + Y{t)dw(t), te[0,T), ( y 

I y(T) = y T . [2A) 

Assume that 1/^(0; M n ) is a separable Hilbert space. For any N G N, write = {ti \ t% = 
JfrT, £ = 0, • • • ,2^}. For any k G {0, • • • , 2^ — 1}, define a sequence of simple processes 

{e«(-, -)},"r by 

= { *.'->W^f >■ (2.2) 
I X[t k ,T]{t)h ki {uj), k = 2 JV - 1, 



where {Mq ; n, M^jv, • • • ,M 2 n_i N } is an increasing sequence of integers. Since Ljr (J);R n ) C 
Lj? (Q;M. n ) (for any < t < s < T), we may assume that {h k {} satisfy the following: 



1) For any fixed k £ {0, • • • ,2^ — 1}, {h ki }^ N is an orthogonal set in L% (0;K n ), and the 



k 



norm \h ki \ L 2 (n Rn) = J^r, and hence \e ki \ L 2 (n L 2 (0 T Rn)) = 1 



^1 



(Q;M") — Y T"' 41111 nelicfc; l e fc«lL|(C>;L 2 (0,T;IR«)) 



2) If < A: < I < 2 N - 1, then {h ki }^ N C {M^f; and 

3) If s = for some iV G N and fc £ {0, ■ ■ ■ ,2^° - 1}, then s = 2 e ~ N »k T/2 e 6 £fy for 

any £ > N . For £ > JV , write k t = 2^ N °k Q . Then, {h kji }™Jj' j C {/i^}£f' for iV < j < £, 

and (J {^fe € i}j =1 ^' is an orthogonal basis of l?-p (f2;IR n ). 
t=N S ° 

Denote by H N the subspace of L§(0; L 2 (0, T; R n )) spanned by {e i}ftf,- ■ • , {^-i^ST -1 '*- 
This is the finite element subspace that we will employ below. Replace L|(0; L 2 (0, T; M n ) by -H^v, 
then the desired numerical scheme follows by trying to find i/n,Yn € Hn such that (jl.2p (with 
/(•,?/(•), Y(-)) replacing by /(•)) holds for rj = and for all u,v € Hpj. 

2 N -lM ktN 

To find the yjv in #iv ; suppose Un = J2 a ki^ki- Choosing u = e k i, v = and 77 = 0, we get 

k=0 i=0 

Zki{t) = f e ki (T)dT, and hence E(z ki (T),y T ) =E/ (zki(r), f(r))dr + Jo ( e ki(r),e tj (T))dT. 

1,3 

Since {efcj} is an orthonormal basis of -ff/v, it follows that E J Q (e^-(r), e k i{r))dT = d k idij- Therefore, 
a fci = prE(/ifci, yr) ~ E y ((r A tfc + i - r A t fc )/» w , /(r))dr. (2.3) 

2^-1 M fejJV 

Similarly, suppose Ijv = £ Pki e ki- By taking u = 0, 77 = and « = e k i to get a corresponding 

k=0 i=0 

-Zfci(i) = /o e ki( T )dw(T), we find that 

/3fci = E((w(tfc +1 ) - w(t k ))h ki ,y T ) - E / ((to(r A £ fe+ i) - u>(t A t k ))h ki , f(r))dr. 

Jo 

We now show the convergence of the sequence {(vn,Yn)} of numerical solutions constructed 
above. 

Theorem 2.1 Let (y,Y) be the transposition solution of \2. Then y^ and Y^ are projections 
of y (viewing as an element of L^iVL; L 2 (0, T; R n ))j andY onto Hjy, respectively. 

00 

By the construction of H]y, it is clear that Hjy C Hm provided N < M. Moreover, Hjy is 

JV=1 

dense in L|(f2; L 2 (0, T; W 1 )). Hence, 

,T f T 

E/ \vn(t) -y(r)\ 2 dT + K / |Yjv(r) - Y(r)\ 2 dr -> 0, as JV -> 00. (2.4) 
Jo JO 



Furthermore, starting from (|2.4j) . we can show the following convergent result. 
Theorem 2.2 As N -> oo, (y N ,Y N ) tends to (y,Y) in L§(ft; Z>([0, T], M n )) x L|(fi; L 2 (0, T; M n )). 
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